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Abstract 



We argue that the non-perturbative Schwinger mechanism may play an im- 
^ , portant role in the hadronic production of charm. We present a flux tube 



model which assumes that the colliding hadrons become color charged be- 



^ , cause of gluon exchange, and that a single non-elementary flux tube is built 

H : 

up as they recede. The strong chromoelectric field inside this tube creates 
quark pairs (including charmed ones) and the ensuing color screening breaks 
the tube into excited hadronic clusters. On their turn these clusters, or 'fire- 
balls', decay statistically into the final hadrons. The model is able to account 
for the soft production of charmed, strange and lighter hadrons within a uni- 
fied framework. 
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I. INTRODUCTION 



Charm production in hadronic and nuclear collisions is presently a subject of considerable 
interest, as charmed particles are expected to be copiously created in relativistic heavy 
ion reactions at BNL's Relativistic Heavy Ion Collider (RHIC) and CERN's Large Hadron 
Collider (LHC). On one side, high energy electrons and muons coming from semileptonic 
decays of these charmed particles will provide a significant 'background' which may obscure 
dilepton signatures of a quark-gluon plasma On the other side, charm production itself 
may be an interesting probe of the QCD plasma due to the mass scale of c quarks. It has been 
suggested that the observation of enhanced charm production in heavy ion collisions could 
provide a way for measuring the temperature of a hot gluon plasma or the thermalization 
time of the initial partonic system 0. 

In order to study charm production in relativistic heavy ion reactions one needs to 
have a good understanding of the production mechanisms operating in nucleon-nucleon 
collisions. Only in this case one can make meaningful extrapolations from hadronic to 
nuclear collisions, and use them to identify new phenomena associated with dense matter 
formation. Charm production in hadronic collisions has been extensively investigated with 
perturbative QCD. Of particular interest for such studies are the recent calculations of 
heavy quark hadroproduction carried out to order by various authors [Q. These next- 
to-leading-order calculations seem to describe very well several experimental features of 
bottom production, but when applied to charm they give a cc cross section which depends 
strongly upon the choice of the renormalization-factorization scale. This renders it difficult 
to assess from direct comparisons with data whether hard processes really represent the 
dominant source of charm in hadronic collisions. Furthermore, there are important features 
of the experimental data which are not reproduced by these perturbative calculations. Some 
examples are the small back-to-back pt correlations found in the charm-pair distribution , 
or the enhancement of leading charm production observed at large 

This leaves open the possibility that non-perturbative mechanisms have an important 
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role in charm production. In fact, the relatively low mass of the c-quark places the hadronic 
production of charm on the border between perturbative and non-perturbative phenom- 
ena. Therefore, approaches from the non-perturbative side may prove instructive. Such 
approaches are usually in the form of phenomenological models. Although less rigorous, 
they provide useful physical images whose parameters, like the string constant or the vac- 
uum pressure, should be object of more fundamental theories. However, the usual mecha- 
nisms for soft particle production do not give a good description of charm hadroproduction. 
The commonly used string models predict negligible charm yields, because the tension of a 
quark- ant iquark flux tube is too small for it to break into a cc pair ||^. Statistical (hydro- 
dynamical) models also do not predict any significant charm production, as typical fireball 
temperatures are much smaller than the charm mass. 

It is then somewhat surprising to find out that by properly blending the string and sta- 
tistical approaches one is able describe very well the hadroproduction of charmed particles. 
Such a hybrid approach to particle production (we call it the "firetube" model) was devel- 
oped in Ref. ^ to study the rapidity and transverse momentum distributions of pions and 
nucleons in hadronic collisions. In the present paper we show how the firetube model can 
be extended to describe charm and strangeness production. 

The general idea of the firetube model is as follows. We assume that the colliding hadrons 
exchange soft gluons when they pass through each other, and as a result acquire color charges. 
The receding hadrons then become connected by a flux tube which confines the strong 
chromoelectric field created by these charges. Quark pairs (including cc) are created by the 
chromoelectric field via the Schwinger mechanism, and the resulting color screening breaks 
the flux tube into lumps of highly excited hadronic matter ( "fireballs" ) which subsequently 
decay thermally into the observed hadrons. The firetube model shares many aspects with 
the usual models based on the fragmentation of a classical string, such as the Lund model 
P]. However, there exist basic differences. First, because of the gluon exchange, in the 
firetube model the string tension can be much larger than that of an elementary string 
between a quark-antiquark pair. Second, the final hadrons come from the thermal decay of 
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fireballs. These points bring some new aspects into the mechanism of hadron production. 
For example, in a common string fragmentation model, the hadrons are produced directly 
from the break up of a quark- ant iquark string via the Schwinger mechanism. Therefore, 
in such a framework the so-called K/ti ratio is closely related to the string constant k. In 
contrast to this, the mechanism for pion production in the firetube model is essentially the 
thermal decay of fireballs, and there is no contradiction here between the small i^/vr ratio 
and a large hadron-hadron string tension. Furthermore, if the string constant is as small as 
the standard value k ~ 1 GeV/fm of a gg-string, the Schwinger mechanism does not produce 
any significant amount of charm, as we have already mentioned. However, the Schwinger 
pair-creation rate is very sensitive to the value of k, and the larger string constants found 
in the firetube model drastically change the picture, providing a mechanism for abundant 
non-perturbative production of charm. 

The aim of this paper is to investigate such mechanism. A brief description of the 
firetube model is given in Sec. and results for pion production are presented. In Sec. ^ 
we discuss charm production in the framework of this model. The total charm cross section 
is calculated for proton-proton collisions and the result compared to experimental data. 
The longitudinal and transverse momentum distributions of charmed mesons and baryons 



are also obtained and compared to measurements. In Sec. [I^ we investigate strangeness 
production and discuss K and A spectra. Finally, Sec. is used for some further comments 
and conclusions. 



II. THE FIRETUBE MODEL 

When two hadrons collide, several sea partons (assumed to be gluons here) may be 
exchanged between them. As a result, these hadrons become colored objects linked by a 
flux tube. Let 5* be the cross section of this tube, which should be of the order of the 
geometrical size of the colliding hadrons. Let also Q be the color charge at the end points of 
the flux tube, measured in units of the elementary color charge Qo of a quark (g^, = JA/3gs, 
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where Qs is the QCD couphng constant). The chromoelectric field E inside the tube can 
be calculated from Gauss' law to be -E = qoQ/S. Placed in such a field, a quark whose 
color charge points into the E direction in color space (we omit SU3 algebra indices for 
simplicity) will experience a force QqE = q^Q/S. Also, the string constant k is related to 
the field energy density by k = E'^S/2. From this we get 

K = KoQ^ (1) 



qoE = 2koQ (2) 

where we have defined Ko = q^/{2S). We identify Ko as the string constant of an "elementary" 
{Q = 1) color triplet string. It is worthwhile to mention that the scaling relations we are 
using — S* ~ ^ E Q^, k, Q"^ — are very different from what one obtains with the 
MIT bag model. In this case the balance between the field energy density and the vacuum 
pressure leads to an increase of the fiux tube cross section with the charge, S* ~ and the 
chromoelectric field and string constant scale as ~ Q^, ~ . However, lattice QCD 
calculations of fiux tubes generated by sources in different representations of the color group 
do predict that the scaling of S, E and k with Q is the one we have used above [^j. 

The end-point charge Q is not necessarily the same for every collision. It fluctuates 
because different numbers of gluons can be exchanged, and also because of the SU3 addition 
of color charges. Assuming that each gluon exchange is a step of a random walk in color 



space [|T0|, it can be shown that the charges generated on the hadrons after the exchange of 
n gluons distribute sharply around the mean value {3/2)y/n, where the 3/2 factor accounts 
for the gluonic octet charge. Thus, k and qoE fluctuate around 

9 

K ^ -KoU , qoE ^ ?>Ko\fn , (3) 

showing that the exchange of even a modest number of gluons gives rise to string constants 
and chromoelectric flelds signiflcantly larger than the elementary ones. 
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In order to determine completely the statistical distribution of the color charge Q we 
must know the probability of having n gluons exchanged in a collision. We assume that this 
is given by a truncated Poisson distribution, 

P„ ~ z/7n! , n > 1 , (4) 

where i/ is a parameter related to the average number of gluons by n = z//(l — e~^). 

The constant chromoelectric field E produces quark-antiquark pairs inside the flux tube 
by a process similar to the Schwinger mechanism of electron-positron creation in QED. 
The qq production rate per unit volume per unit time can be calculated (in the Abelian 
approximation) from Schwinger's formula 

where ruq is the quark mass. Corrections for the final state interaction of the qq pair can 
be introduced into this formula by taking qoE = 2koQ — i^o- In our calculations we also 
corrected Eq. for the finite transverse size of the flux tube, using the semiclassical formula 
of Ref. (see their eq.27). For the cases of interest to us, the simple semiclassical result 
is not very different from the more precise correction of Ref. [^]. In principle we should 
also correct the Schwinger formula for the finite length of the flux tube [T^. But it has 
been shown that if the tube end-points are moving rapidly (which is our case) such 
corrections are suppressed by relativistic dilation effects. For this reason we have neglected 
longitudinal size corrections. 

The quarks created by the chromoelectric field have their color charges aligned with E 
in SU^ space (the other orientations will be ignored as they have a much smaller production 



rate |J0|,|T5|) and tend to screen the end-point charges Q, providing a mechanism for the flux 
tube fragmentation. Although other processes such as collective instabilities of the vacuum 
can be invoked to explain the break up of the tube into fireballs, we will assume that pair 
creation gives the dominant mechanism. This has the advantage of allowing for a simple 
estimate of the firetube fragmentation rate per unit length per unit time as 



a;^|E^. (6) 

where the sum extends over the quark flavors {u, d and s in practice). This is, admittedly, 
a rough treatment of screening effects, but we did check that the typical break-up time we 
obtain in this way is comparable to the collapse time of the color field calculated with quark 
transport models . 

The firetube model has a relatively small number of parameters. For proton-proton 
collisions we take n = 2.0, Kq = I GeV/fm and S = 1.5 fm^. Note that these values are quite 
reasonable: is the usual qq string constant |]8|,p!5|, and the firetube radius R = \J S/n ~ 
0.7 fm is almost the same as that of the proton. A few other parameters define the minimum 
fireball mass (1 GeV), regulate the behavior of the leading particles, and determine how the 
effective temperature and longitudinal expansion rate of a fireball depend upon its mass [0 . 
For the quark masses we use the constituent values m„ = = 300 MeV, = 450 MeV 



and iric = 1.5 GeV 15 



Having defined all this, we can perform Monte Carlo simulations of the firetube formation, 
its fragmentation into fireballs, and the thermal decay of these into the observed hadrons 
0. A typical result of such calculations is presented in Fig.l, where we show the rapidity 
distribution of charged particles (mostly pions) produced in proton-proton collisions at i/i = 
20 GeV and 53 GeV. We see that the model calculation is in good agreement with the 



experimental data [|T6| at both energies. Results of similar quality are obtained for the 
transverse momentum distributions. 



III. CHARM PRODUCTION IN THE FIRETUBE MODEL 

As it can be seen from Schwinger's formula, Eq.(|^), the production rate of charmed 
quarks is extremely sensitive to the value of qoE. For example, if qoE were of the order 
of 1 GeV/fm as usually quoted in the string model, the dominant exponential factor in 
Eq.(^ would be of the order of 10^^^, and no reasonable interaction volume and time scale 
for hadron-hadron collision would account for the observed charmed particle production 
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cross section. On the other hand, in our model the average value of qoE is approximately 
3 GeV/fm and the average value of the dominant exponential factor for charm becomes as 
large as 10~^, for the parameter values given above. 

In order to obtain the charm production cross section we proceed as follows. First, we 
calculate the available space-time volume VT as the total area swept by the flux tube in 
the space-time plane times the tube cross section S. This is completely determined if the 
fragmentation scheme of a firetube is specified. The total number of charmed quarks is 
then twice the value VT x TZc- In Fig. 2 we show the total charm production cross section 
calculated in this way, and compare it to experimental data |lT7|-pO|]. The agreement is seen 



to be good, specially if we have in mind that the same parameter set is used at all energies. 
In our model, the energy dependence of the charm production rate comes essentially from 
that of the space-time volume VT swept by the firetube, which increases asymptotically as 
ln(s). 

To calculate the momentum distribution of charmed particles, we assume that the 
c quarks produced by the mean color field are distributed among the tube fragments with 
probabilities which are proportional to the mass of each fireball. The charmed particles are 
then emitted from the fireballs following their longitudinal expansion and thermal decay, in 
the same way as the other mesons [0]. In Fig. 3 we show our calculation for the xp distribu- 
tion of D mesons produced in proton-proton collisions at ^/s = 27 GeV (solid line) together 
with experimental points The agreement is very satisfactory. The transverse momen- 
tum distribution of charmed mesons predicted by the model is also in a good agreement 
with experimental data [T^, as shown in Fig. 4. 

In addition to the particles coming from fireballs, hadronic spectra also get a contribu- 
tion from the leading particles. The mechanism for tube breaking into fireballs assumes 
a fragmentation rate u homogeneous in space-time. There is no apriori reason for this to 
remain valid at the two end-points which contain the valence quarks, because of the different 
boundary conditions. In fact, to reproduce the observed leading nucleon spectra we have to 
require that the two proton-like extremities detach from the firetube with a probability rate 
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that is constant on their world-hnes. When these hadrons separate from the firetube, there 
exists a chance that they turn into charmed particles, as the detaching mechanism should 
also be related to the production of qq pairs. Thus, the probability of having a charmed 
leading particle can be estimated as IZd J^qT^g- We have calculated in this way the xp spec- 
trum of leading charmed baryons (A^'s) produced in proton-proton collisions, and the result 
is indicated by the dotted line in Fig. 3. Note that these calculations do predict a substantial 
leading particle effect, as the Ac spectrum is much harder than the D/D distribution. For 
pion-nucleon reactions, where detailed observations of the leading charm effect have been 



performed [||, our calculations reproduce very well the experimental results |2T|. We should 
mention that other non-perturbative effects, like intrinsic charm or valons [Q, have 
also been invoked to explain the leading particle behavior. 

Another important feature of the firetube model is that we expect no strong transverse 
momentum correlations in the D — D pair spectra, as the charmed particles are emitted from 
fireballs which decay statistically. This seems to be consistent with most of the available 
data, which show almost no back-to-back angular correlations in the transverse plane 

A final point concerns charm production in hadron-nucleus collisions. Some preliminary 
extensions of our model indicate that the forward (x^ > 0) charm production cross section 
is essentially proportional to the target mass number A. However, for negative xp the charm 
production rate seems to increases more rapidly than linear in A. 



IV. STRANGENESS PRODUCTION 

We have already seen that the large string tension arising in hadronic collisions can 
accommodate both the pion and charmed meson spectra quite satisfactorily. At this point 
one may worry about strangeness production. For simplicity we will assume here that 
strange particles are created solely through the Schwinger mechanism, in the same way as 
charmed particles. This means we are ignoring thermal production of strangeness, which in 
principle could take place inside the fireballs (for charm this process is certainly negligible 
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because of the large mass of the c-quark). In Fig. 5 we show the xp spectrum of neutral 
kaons produced in proton-proton collisions, calculated exactly as we did for the D mesons. 
We see that the agreement of our results (shown as a solid line) with the experimental data 



2^ is very reasonable. Transverse momentum distributions predicted by the model are also 



in good agreement with experimental results. The production of strange leading particles 
can also be calculated along the same lines followed for the charmed ones. Our result for the 
spectrum of leading A baryons produced in p-p reactions is shown in Fig.5 (dashed line). 
The accordance with experimental data is again quite good. 



V. FINAL REMARKS 

In this paper we have discussed the possibility that a non-elementary flux tube is created 
between two colliding hadrons because of gluon exchange. This picture is consistent with 
charmed particles being produced at large rates by the color field inside the tube, via the 
Schwinger mechanism. We developed a simple model based on this idea and showed that it 
could account quantitatively for the production of charm in proton-proton reactions. The 
same model also describes very well the creation of lighter hadrons such as pions and kaons. 
It is a noteworthy aspect of our results that a successful description of the pt and xp spectra 
of such different particles as vr's, K^s and D's can be achieved in a unified manner. 

An important issue refers to the competition between the mechanism for non-perturbative 
charm production proposed here and perturbative partonic processes. Hard processes cer- 
tainly contribute to the charm yield but, as mentioned in the introduction, the present 
uncertainties in perturbative calculations make it difficult to determine the size of such a 
contribution. Even though we have shown that non-perturbative production by the color 
mean field can easily account for the observed charm cross section, we do not claim here 
that this mechanism represents the only significant source of charmed particles. It should be 
clear that the precise values of our charm cross section also suffer from uncertainties coming 
from parameter values. For example, by adopting a slightly larger value for the mass of 
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the c-quark we would reduce our charm production cross section without changing any of 
the model predictions for non-charmed particles. The basic point we want to stress here 
is that, within a reasonable set of parameters, our model shows that an important part of 
the charmed particles produced in a hadronic collision may have a non-perturbative origin. 
So, from our point of view, soft and hard mechanisms may well be both necessary in order 
to account for the bulk of charm production. What is certain is that several experimental 
observations cannot be described on purely perturbative grounds, and require a signifi- 
cant contribution of soft mechanisms in order to be understood. In this paper we presented 
a consistent non-perturbative framework that can describe not only these but most features 
of charm production. For the reasons discussed in the introduction, it is of interest to inves- 
tigate how the mean field mechanism for charm production discussed here scales to heavy 
ion collisions. Work along these lines is in progress. 
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FIGURES 

FIG. 1. Rapidity distribution of charged particles produced in p-p collisions at ^/s = 20 GeV 



and 53 GeV. Experimental data are from Ref. [16|. The curves represent the result of Monte-Carlo 
calculations with the firetube model. 

FIG. 2. Total cross section for charm production in p-p collisions as a function of cm. energy. 
The line is the result of the firetube model. Experimental points (solid squares) are from Refs. 
1 17, 18 1 . The open squares refer to data obtained with proton-nuclei reactions ]l9| , pm . 



FIG. 3. Feynman xp distribution of charmed particles created in p-p reactions at Pj„c = 
400 GeV/c. The solid and dashed curves correspond to our calculations for D/D mesons and Ac 
baryons, respectively. Data points for D/D are from Ref. [17|. 



FIG. 4. Transverse momentum distribution of charmed mesons for the same reaction as in 
Fig.3. 

FIG. 5. Feynman xp distribution of strange particles produced in p-p collisions at Pine = 
205 GeV/c. The solid and dashed curves correspond to our calculations for the K^/K^ mesons 
and A baryons, respectively. Experimental data are from Ref. 124]. 
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